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ABSTRACT 

Spectra of six WN + OB Wolf-Rayet systems obtained with the IUE 
are analyzed for phase-dependent variations. Periodic variability 
at emission-line frequencies is detected in V444 Cyg, HD 90657, HD 
211853, HD 186943 and HD 94546 on low dispersion SWP (XX 1200 - 1900 A) 
images. No changes in the low dispersion spectra of HD 193077 are 
apparent . 

We find the variations in the UV to be similar in nature to 
those observed in optical spectra of various WR sources. That is, 
there is a strengthening of absorption components in P Cygni - type 
features at orbital phases in which the 0-sta,r is behind the WR wind. 
With the aid of a computer code which models this fvpe of variations, 

and through a comparison with HD 193077, the dominant mechanism 

producing the variations is shown to be selective atmospheric 

eclipses of the 0-star by the WR wind. Based on this interpretation, 

a straightforward technique is applied to the line of N IV 1718, 

by which an optical depth distribution in the WN winds of the. form 

T *r"l is derived for 16 < r < 66 R . 

a 

The presence of the N IV ionization stage out to at least 
66 is consistent with the proposed constant ionization structure 
in WN winds (Willis, 1982). 

Phase-dependent variations in the width of the C IV 1550 absorp- 
tion component in V444 Cyg, HD 90657 and HD 211853 are interpretated 


as wind-wind collision effects. 



TABLE OF CONTENTS 


Page 

ABSTRACT iii 

LIST OF TABLES v 

LIST OF FIGURES . vi 

ACKNOWLEDGEMENTS 

Chapter 

I. . INTRODUCTION 1 

II. IUE OBSERVATIONS AND INTERPRETATION 9 

Introduction 9 

Observations • ••••• 10 

Data Reduction * 13 

Ultraviolet Variations 23 

Interpretation 29 

Wind-Wind Collision 41 

Possible Fe V and Fe VI Contributions. ... 46 
Radial Velocity Curves ..... 49 

III. EMPIRICAL OPTICAL DEPTH DISTRIBUTIONS 52 

IV. CONCLUSIONS AND RECOMMENDATIONS. 69 

REFERENCES * 72 

APPENDIX: MODEL PARAMETERS . 76 


LIST OF TABLES 


Table Page 

1 Observed Systems IX 

2 Observed Targets and Measurements . . , 16 

/ 

3 Averages and Standard Deviations 22 

4 Suggested Line Identification 48 

5 Impact Parameters and Computed Optical Depths 56 




v vi 

LIST OF FIGURES 

Figure Page 

1 IUE Spectrum of HD 193077 ... 14 

2 IUE Spectra of 

a) V4444 Cyg 24 

b) HD 211833 25 

c) HD 90657 26 

d) HD 186943 27 

e) HD 94546 28 

3 Geometrical Setting . 31 

4 Equivalent Widths of N IV 1718. . . * ,36 

5 Comparison of Model and Observed Variations 38 

6 Ratios of Spectra at Orbital Phases "0.0 and "0.5 

a) V444 Cyg, HD 211853 and HD 90657 42 

b) HD- 186943 and HD 94546 / 3 

7 Radial Velocity curves at N IV 1718 50 

8 Optical Depth Distribution 55 

9 Possible Changing Excitation Structure 67 

10 Grid of Model (a,e) Values. 80 


ACKNOWLEDGEMENTS 


This thesis represents the final stage of my formal academic 
instruction. Many persons have played important roles throughout 
this formative period by introducing concepts, influencing 
philosophical perceptions, and providing valuable instruction. Their 
guidance, interest and support have made this thesis possible. 

In a chronological order, I must first thank Franz. Kaenigsberger 
who sparked in me the curiosity and interest in celestrial bodies 
and in natural phenomena, and Linnea Horowitz who set the stage for 
my first years of formal education. To both my parents, I express 
my graditude. 

I am deeply indebted to Anthoni Aveni, whose advice and 
encouragement may be the single most important influence responsible 
for my career in Astronomy. I sincerely thank him for his continuing 
interest in my progress during the past eighteen years. 

I am profoundly grateful to Manuel Peimbert for the valuable 
guidances and support he has always given me. 

Claudio Firmani introduced me to astronomical research, and 
is responsible for the framework within which most of ray research has 
been set. I gratefully acknowledge his influence on my work. 

Roberto Avena's persistence, confidence and backing are largely 
responsible for my graduate degree. His support has also enabled me 


viii 


to gain important experience at major centers of research. I am 
fortunate to have his collaboration in all my endeavors. 

Lawrence Auer hat- been more that a thesis adviser, he has been 
a true friend, generous with his time, and willing to share his vast 
experience and knowledge. His endurance of lengthy long-distance 
telephone discussions, and his patient and unrelenting efforts to 
convey to me an intuitive understanding of stellar atmospheres are 
truly appreciated and gratefully acknowledged. 

I thank Phil Massey and Anthony Moffat for useful discussions 
and for reading and answering my lengthy letters. Many of the 
arguments used in this thesis originated in this correspondance. 

My sincerest thanks to David van Blerkom for useful discussions 
and for his role in laying the foundations for this investigation. 

This thesis would not have been possible without the IUE 
observations on which it is based. We are grateful to the IUE 
staff. Skip Schiffer, and especially Robert Panek, who provided 
invaluable assistance throughout the data acquisition, reduction 
and analysis phases. Thanks are also due to Wayne Warren and the 
National Space Science Data Center far the IUE spectra drawn from 
their archives . 



I thank Virpi Niemela for providing results of her investi- 
gation prior to publication, and Jose Franco for helpful discussions. 


ix 


Invaluable research experience was acquired while working 
with Jean Swank at the Laboratory for High Energy Astrophysics, 
NASA/GSFC. I am grateful for her support and interest. 

I thank Gordon Garmire for useful comments on an earlier draft 
of this work, and for editing the final version. I also thank Satoshi 
Matsushima, Don Page, Lawrence Ramsey and Daniel W'eedman for useful 
suggestions. The corrections and the production of the final 
version were tremendously simplified thanks to AJ3E, created by 
Lawrence Auer. 

Roberto Avena and Alberto Garcia are thanked for the drawings, 
as is Ruth Hollinger for her assistance in dealing with official 
matters. 

Financial support for graduate work was provided by the 
National Science and Technology Coqncil of Mexico (CONACYT) and 
is gratefully acknowledged. The IUE observations and part of the 
data reduction were supported by a grant from NASA. The. Instituto 
de Astronomia of the Universidad Nacipnal Autonoma de Mexico 
provided office space and computer time during the final stages of 
this investigation. 




1 


Chapter X 
INTRODUCTION 

Since their discovery in 1867 by C. J. E. Wolf and G. Rayet, 
the Wolf-Rayet (WR) stars have remained as one of the most interesting 
groups of stellar objects. Although considerable theoretical and 
observational work has been dedicated to them, and although advances 
have been made in determining their properties, an understanding of 
these objects has not been achieved. 

For the more recent descriptions of their properties and 
theoretical advances the reader is referred to the four symposia 
and colloquia dedicated to WR stars and related objects (Gebbie and 
Thomas, 1968; Bappu and Sahade 1973; Conti and de Loore 1980; 
de Loore and Willis 1982) , and references cited therein. We will 
summarize only the characteristics of WR stars needed to set the 
framework within which the present investigation is contained. 

The spectra of WR stars are characterized by broad and 
intense emission lines. On the basis of these emissions lines, 
two general categories or subtypes have been established: The WN's 

which show mainly lines of nitrogen and the WC's where carbon and 
oxygen lines prevail. In both cases, helium lines are present, 
and hydrogen lines are weak or absent. Both categories are 
subdivided according to degree of ionization. The WN's are 
classified according to decreasing degree of ionization from WN2 
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to WN10. Analogously, WC’s ate classified from WC5 to WC9. In 
most cases, for a given star, a large range in ionisation and 
excitation conditions are represented, the most extreme case 
having four stages of ionization of the same element visible in 
its spectrum (Massey and Conti 1980) . 

Many of the emission lines exhibit P Cygnl type profiles. That 

is, the emission is accompanied by a shortward-shifted absorption 

component, which is characteristic of envelopes or winds in 

expansion with respect to a continuum source. It is the presence 

of these absorptions and the large widths of the emission lines 

that prompted Beals (1939) to first suggest the existence of 

massive outflows of material from WR stars. The emission line widths 

are generally of the order of 3000 ,km s”^, indicating outflow 

—1 

velocities of at least 1500 km s . 

Two fundamental problems concerning WR stars are stiXL 
unsolved: Their evolutionary status and the mechanisms powering 

the enormeous mass-loss rates observed. Although WR stars are 
generally considered to be in the late stages of the evolution 
of massive stars (i.e., post Main Sequence), the conditions necessary 
for the progenitor to evolve into a WR are unclear. Abundance 
analyses indicate that WR's have nuclear processed material at 
their surfaces, with WN subtypes corresponding to stars in which 
the products of the CNO cycle are on the surface (and hence in 
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the wind) and WC's showing products of helium burning (Willis 
and Wilson 1978; Paczynski 1973). Only for the most massive 
stars are mass-loss rates through stellar wind sufficient to 
strip away enough layers to uncover nuclear processed material. 

In the remaining cases, mass-loss through Roche Lobe overflow in 
binary systems, extensive mass-loss during the red supergiant stages, 
or some other means of mass-loss is necessary. 

Considerable work has been dedicated to the construction 
of evolutionary models leading to WR stars. These Include the 
effects of mass-loss during the Main Sequence, binary membership, 
rotation, Roche Lobe overflow, and combinations of these. The 
reader is referred to the numerous publications cited in the 
reviews by Chiosi (1982) and de Loore (1982) . 

In addition to the mass-loss necessary to generate a, WR star 
from a massive progenitor, the WR's are currently in a phase of 
moss-loss via a stellar wind. The mass-loss rates which have been 

• _ e _ 1 

deduced are M > 10"“ yr , and are apparently independent of 
subtypes and stellar parameters (Barlow, Smith and Willis 1981) . 

The mechanisms responsible for this large rate are not fully 
understood. Apparently, radiation pressure can only account for 
part of the mass outflow (Abbott 1982). Thus, other mechanisms 
have been invoked in which mechanical energy is generated and then 
transmitted to the wind, powering the flow. For example, the 
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presence of a hot and thin corona just above the photosphere has 
been proposed (Thomas 1973} Hearn 1975; Cassinelli, Olson and 
Stalio 1978). Other models include a two-component gas, in 
which a population of radiatively-driven condensations plough 
through and are confined by ram pressure of, an ambient non- 
radiatively driven r * (Lucy and White 1981; Lucy 1982). Finally, 
Sreenivasan and Wilson (1982) have proposed a model in which a WR 
star can be viewed as a rapidly rotating core which is being slowed 
down by the extended envelope with which it is not in co-rotation. 
The shear energy is then transported to the wind and provides 
mechanical heating. 

r ‘.UMy ly, an understanding of the structure of WR winds can 

preside an insight into the driving mechanisms. By wind structure, 

one means the functional dependence of the velocity, density and 

temperature on position throughout the wind. Due to conservation 

of mass flow, the density and velocity at each point in the wind 

2 

may be related (ft * r v(r)p(r)), and thus if the velocity 
distribution can be determined, the density distribution follows, 
and vice-versa. 

Several attempts have, been made to derive the structure in 
WR winds by modeling the infrared continous energy distribution 
assuming free-free emission (Hartmann and Cassinelli 1977; van der 
Hucht 1979; Hartmann 1978). These resulted in density distributions 
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corresponding to a wind with an extended constant velocity 
distribution. Consistent with these results is the analysis of 
continuum eclipses of V444 Cyg, a well-known WN + OB binary system, 
where a constant velocity flow is Inferred for radial distances 
r > 12 R (Cherepashchuk, Eaton and Khaliullin 1981). If this 

O 

result is correct , then the wind must undergo the initial 
acceleration close to the continuum-emitting core, but there is 
no information regarding this region. 

Perhaps the only means of determining wind structure is through 
the analysis of emission line profiles. The basic picture of a WR 
star consists of a hot continuum-emitting stellar core surrounded 
by a spherical, rapidly expanding extended atmosphere, in which 
the emission spectrum is assumed to originate. 

The treatment of the radiative transfer in moving 
atmospheres has been given by Sobolev (19 53, 1960) and extended 
by Castor (1970a) . Sobolev was the first to recognize that the 
presence of a velocity gradient in the extended envelopes of stars 
such as WR's actually simplifies the line transfer problems 
for it dominates the photon escape and thermalization process, and 
implies a geometric localization of the source function. That is, 
in a rapidly moving extended atmosphere with an expansion velocity 
far in excess of the thermal velocity, radiative interaction in 
the lines with distant parts of the atmosphere is negligible 
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and thus the line transfer is locally constrained. Hence, a line 
photon emitted in a certain direction at a given frequency in a certain 
place in the wind is either absorbed immediately, or escapes, and its 
frequency corresponds to the line-of-sight velocity of the emitting 
material. The geometrical region from which the emission at any one 
frequency arises is a surface of constant radial velocity. These 
surfaces extend over large regions , and their shape depends upon the 
nature of the velocity field. Thus, we can not associate a given 
frequency in the line profile with a specific position in the wind, 
but only with a line-of-sight velocity. Insofar as large regions 
of the wind may have the same line-of-sight velocity, the wind 
structure cannot be extracted directly from the observed profiles. 

The technique usually employed Is to model the profiles, using the 
Sobolev Theory, and to apply different source functions and optical 
depth distributions. However, this leads to ambiguous results, 
since line, profiles can be reproduced with a variety of dissimilar 
combinations of optical depth and source functions (Castor 1970 a) > 

Thus, as stated by Mihalas (1978, p. 568), a thorough analysis of 
the spectroscopic data, with the goal of diagnosing physical conditions 
in the flow semiempirlcally is necessary and might prove to be very 
rewarding. 

Up to now, all the attempts to model profiles in WE stars 
have been concentrated to fitting profiles observed in single stars; 
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i.e., those not associated with o* companions. It has long been 
known (c. f . Wilson 1940) that, line profiles in WR binary systems 
undergo periodic phase-dependent variations. Hence, the interaction 
effects suggested by these variations .represent an added complication 
to an already difficult problem, especially since the mechanism 
responsible for the variations are not fully understood. However, 
binary systems can provide an extremely valuable quantity: the 

scale length within the wind. That is, the orbital separation, 
which is presumably known in the well studied systems, can be used 
as a yardstick to establish the geometrical location of emitting 
and absorbing material within the wind. Specifically, thirty- three 
years ago Munch (1S50) proposed that the eclipses of the OB component 
by the WR wind are largely responsible for variable P Gygni-type 
profiles observed in these binary systems. Recent observations at 
ultraviolet (Willis et al. 1979; Eaton, Cherepashchuk, and Khaliullin 
1982) as well as optical (Cowley et al. 1971; Massey and Niemela 
1980) frequencies seem to indicate that selective atmospheric 
eclipses are indeed the dominant mechanism. If this is the case, then 
the OB component in WR binary systems can be used as a probe of the 
wind opacity distribution, and hence restrict the number of free 
parameters in model profile computations. 

In this thesis we will investigate the possibility of using 
phase-dependent profile variations in WR + OB binaries in order to 
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derive the wind structure of WR stars. We will show that the major 
profile variability results from selective atmospheric eclipses, and 
we will empirically derive a first order optical depth distribution for 
the WR winds. 


9 


Chapter XI 

IUE OBSERVATIONS AND INTERPRETATIONS 
Introduction 

The ultraviolet (UV) spectral region of WR stars has been 
extensively observed since the advent of UV astronomy. References 
to the early results are listed by Johnson (1978) , while subsequent 
observations are reported by Willis and Wilson (1978) and van der 
Hucht (1978). However, until the launch of the International 
Ultraviolet Explorer (IUE), UV observations of WR stars were 
limited to the brightest objects^, and to low resolution spectroscopy. 

The IUE was launched in early 1978, becoming the first astronomical 
satellite to be placed in a geosynchronous orbit, permitting a 
greater accessibility of targets. It carries two UV echelle 

o 

spectrographs covering the wavelength region 1150 - 1950 A and 
0 

1900 - 3200 A, respectively, permitting the acquisition of high 

O 

resolution (0.1 - 0.2 A) UV spectra for the first time. Low 
resolution (-6 A) spectra can also be obtained of much fainter 
objects than previously possible. The description and performance 
of the IUE have been summarized by Boggess (1978, 1979). 

The first year of IUE observations of WR stars has been 
summarized by Nussbaumer et al. (1979) , and subsequent observations 
are reported in two volumes published by NASA (Chapman 1981; Kondo, 

Mead and Chapman 1982) . 
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The data obtained from ground-based observations suffer from 
a fundamental limitation: all of the lines observed are subordinate 

transitions that arise from levels with high excitation potentials, 
and hence with small populations outside the regions of high 
temperature and density. The column density of absorbers in these 
lines is therefore quite small, and one observes only the innermost 
layers of the wind, just outside the photosphere. In contrast, 
the UV spectral region contains resonance lines arising from the 
ground states of the dominant ionization stages of abundant elements. 
The column densities in these lines are so great that one may 
sample the outermost parts, of the wind. In addition to the 

t 

resonance lines, the UV region also contains subordinate transitions 
that arise from levels of low excii ,tion potentials. 

In this chapter we will report the results of our IUE 
observations of six WN + OB systems. 

Observations 

Table 1 lists the Wolf-Rayet stars selected for this study, 
along with the available orbital parameters. The systems all contain 
an early WN component (WN4-6) , have short orbital periods, are 
bright enough for short exposure times, and most are known to undergo 
periodic profile variations in their optical spectra. Two of the 
systems, V444 Cyg and HD 211853, are eclipsing binaries, though the 
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latter is apparently part of a quadruple system (Massey 1981) . 

Thus, orbital elements for this system may be uncertain. HD 193077 
is a proposed triple system, containing a WN6 + compact object in a 
close orbit (P * 2.3 days), both of which are in orbit about an OB 
star at a considerable distance (Lamontagne et al. 1982). 

The observations were carried out with the IUE at NASA/GSFC with 

Q 

the Short Wavelength Prime (SWP) camera (XX 1200 - 1900 A) , through 
the large aperture and in tne low dispersion mode. Exposure times 
were chosen so as to optimize the count levels in the continuum and 
in the N IV 1718 emission line. Hence, He II 1640 was just satu- 
rating on most exposures. Low dispersion images were decided 
upon in order to maximize phase coverage of the six systems, reduce 
the data reduction complications, obtain good signal to noise 
images of the weaker targets, and permit accurate positioning of 
the continuum. In addition, lines in WR stars are very broad 
(~300Q Km s~l) and are thus ^3 resolution elements wide. 

The six targets listed in Table 1 were observed on each of 
six observing shifts, and an additional 21 SWP images of V444 Cyg 
were drawn from the archives of the National Space Science Data 


Center. 
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Data Reduction 

Data reduction was conducted at the Regional Data Analysis 
Facility (RADF) , NASA/GSFC. Standard IUE procedures were employed 
to obtain absolute fluxes, equivalent widths and wavelength 
positions. Additional flux and wavelength measurements were made 
directly from calcomp plots. Recorded exposure times of trailed 
images were divided by 1.07 to correct for the adjusted trail 
rate (Panek 1980) . No correction for interstellar reddening has been 
applied. 

Since we intend to analyze the variability, it is important 
to discern fluctuations due to -the detection and reduction pro- 
cedures from the actual spectral variations. In addition, we require 
an estimate of the uncertainties in the measurements of velocity, 
flux and equivalent width. 

As pointed out by Ayres (1982) , a sequence of several images 
exposed to the same count levels can be used to determine the standard 
deviation of the individual fluxes about the mean, thus taking 
automatically into account errors resulting from the original image 
processing. We will apply this procedure to derive a quantitative 
estimate of the uncertainties in the measurements, and we will use 
the eight images of HD 193077 to do so. Given our phase, coverage, 
practically no variability in this source was detected in our low 
dispersion spectra. Thus, it will serve as a ''standard” for com- 
parison with the other more variable sources. 
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FIGURE 1 


IUE spectrum of HD 193077 non-dereddened„flux is 
is given in units of 10“^ erg cm" 2 s" 1 * A"^-, and 
wavelengths in Angstroms. 
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In Figure 1 we show the spectrum of HD 193077, obtained 

by co-adding the eight images listed in Table 2. The most prominent 

emission features are marked. The large concentration of emission 

0 

lines is evident in the region X < 1550 A. Clearly, any analysis 

of wind structure requires an unblended line, one which arises in 
the WR wind, and sufficiently intense for reliable measurements 
to be made. The most promising feature in this respect is the N IV 
line at 1718 4* 

In Table 2 we list all the measured quantities for each of the 

IUE images, identified by their SWP (’’Short Wavelength Prime" camera) 

numbers, listed in column 1. SWP numbers smaller than 15500 

correspond to images drawn from archives cf the National Space 

Science Data Center. An asterisk indicates a trailed image. In 

column 2 we list the orbital phase computed from the ephemerides 

in Table 1. Columns 3-5, 8 and 9 contain the positions (in 

units of 10 2 km s“^) of the wings and the extrema in flux of the N IV 

1718 P Cygni feature, with the following notations: V^ and V^ 

correspond to the blue and red wings, respectively, of the absorption 

component: V^ to the red wing of the emission component; and V^ and 

V E are the minimum in the absorption and the maximum in the emission 

respectively. Columns 6 and 10 contain the eye- interpolated mono- 

chromatic flux (in units of 10“ erg cm“^ s”^ A -1 ) of the absorption 

(F.) and emission (F_) extrema, respectively. Columns 7 and 11 
A E 


TABLE 2: Observed Target, and Hea.ure.ent. 
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TABLE 2; Observed Targets and Measurements (continued) 
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|5645* 0.26 -26. ~}l. - 2. 4*. 2 2.0 + 7. 27. 20.4 3.5 


TABLE 2; Observed Targets end Haaaureoents (continued) 
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Notes to Table 2 


Asterisk Indicates trailed exposure. 

Phases computed from periods and initial epochs in Table 1, the 
uncertainties in these phases are: 0.002 (7444 Cyg) , 0.02 (HD 

90657), 0.001 (HD 186943), 0.003 (HD 211853), >0.3 (HD 94546). 
Velocities in units 10^ km s“^. 

Fluxes in units of 10"^ erg cnf^ 3 “1 A”*. 

O 

Equivalent widths in units of Angstroms. 

Phases in parentheses have large uncertainties. 
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contain the corresponding equivalent widths, and finally, column 
12 lists the monochromatic flux in the continuum region near the N IV 
1718 line, obtained by interpolating the surrounding continuum by 
eye. 

In the 38 images obtained during our observing run, we also 

0 

measured flux in the continuum region around 1825 A, which is 
relatively free of lines. The averages + standard deviations for this 

quantity ( F 182 5 ^ » V A> V E» F A» F E» W A» 311(1 W E are listed ^ Table 3. 

The sign indicates percent variations; i.e., (3.d.)/mean. 

According to the results in Table 3, we can state that, with 
respect to HD 193077, there is significant line flux and equivalent 
width variability in V444 Cyg, HD 90657, HD 211853, and moderate 
variability in HD 188943 and HD 94546. On the other hand, little or 
no continuum variability is detected. It is interesting to note 
the large difference between the velocity standard deviations in 
HD 193077 and those of the reset of the sources. Specifically 
we find the standard deviation of relative measurements to be 
~90 km s' 1 from the HD 193077 values. This illustrates the remarkable 
stability of the IUE wavelength calibration procedures, and quantifies 
the uncertainty in establishing the position of extrema in the P 
Cygni profiles. Thus, although velocity measurements are limited 

o 

in general by the "6 A resolution of the low dispersion spectra, 
position of extrema (i.e., line centers) can be determined with much 


\ 
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TABLE 3; Averages and Standard Deviations 
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higher precision. Furthermore, the stability of the wavelength 
calibration permits relative flux measurements at a given 
wavelength to be obtained with considerable reliability. We 
emphasize that our present analysis will require only relative 
flux and wavelength measurements. The velocity standard deviations 
in all the other sources are -200 - 300 km s”^, comparable with 
typical orbital velocities reported, for WR stars. 

Ultraviolet Variations 

Figures 2 (a) - (e) illustrate the UV spectra at individual 
orbital phases of the five variable sources. The variations can 
be summarized as follows: 

o 

a) The P Cygni structure of N IV 1718 A undergoes a 
transition from a strong emission accompanied by a weak 
or absent absorption component when the WR is in front 
of the O-star (<}> - 0) to weak or absent emission with a 
shortward -shifted , very intense absorption component 
at the opposite phase. 

b) He II 1640 develops a shortward- shifted absorption 
component at <p ~ 0, but no absorption is apparent at 
other phases. 

o 

c) The apparent continuum region shortward of -1500 A is 
significantly depressed at <f> - 0. In the more extreme 


FIGURE 2(a): IUE spectra of V444 Cyg. Spectra at similar orbital phases 

have been co-added, and are identified by these phases. 





25 


cmsQfim facs a 

OF POOR QUALITY 



FIGURE 2(b): Individual IUE spectra of HD 186943 Identified by orbital 

phase. 
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FIGURE 2(d); Individual tUE spectra of HD 211853 identified by orbital 
phase. 
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1360 . 1523 . 1600 ., 1840 . 2003 . 


FIGURE 2(e): Individual IUE spectra of HD 94546 identified by orbital 

phases computed with data in Table 1. By analogy with 
the previous systems, we have placed the spectrum with 
strongest absorption (labeled with phase 0.43) on the 
bottom 
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cases this constitutes a variation of 45 Z in flux 
level, while no variations greater than 8% occur in 
the 1825 continuum region in the same Images. 

d) The emission component of C XV 1550 shows a decrease in 
intensity at <J> -0, but at the same time the absorption 
component is narrower instead of broader and deeper, as 
occurs with the N IV 1718 line. 

As already mentioned, the strongest variations are seen to 
occur in V444 C.yg t , HD 90657 and HD 211853, though a more complete 
phase coverage of the^other two sources may disclose similar degrees 
of variability. There is no evidence for continuum eclipses in any 
of the sources except V444 Cyg,' where the additional images obtained 
at $ * 0 display variations consistent with the 30% ecplises observed 

o 

in. the optical region (XX 4200 - 7500 A; Cherepashchuk and Khaliullin 
1973). The continuum eclipses are significant at <P * 0.0 + 0.03. 
Thus, the variations at emission line frequencies occur over a wider 
phase interval than the continuum eclipses, as would be expected from 
selective atmospheric eclipses. 

Interpretation 

Basically, we can describe the observed UV variations as 
the presence of strong absorptions at <f> - 0 and their decrease 
or disappearance at the opposite phase. This is exactly the same 


? 
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behavior that has ben reported for the optical region in a few 
WR binaries (c, f. Munch 1950; Cowley et al. 1971), and we will 
establish a suitable explanation for these variations in terms 
of selective atmospheric eclipses. 

According to this interpretation, as the 0-star passes ' 
behind some portion of the WR wind, it3 continuum radiation is 
absorbed by ions in the wind. The absorption occurs at frequencies 
corresponding to atomic transitions with favorable transition 
probabilities whose lower levels are adequately populated. 

In this respect, these frequencies are selected for absorption, while 
others are not affected. However, because of the large expansion 

i 

velocities in the winds, a relatively wide range of frequencies will 
be affected, as we will explain below. 

For a given transition, each frequency that we observe in an 
emission line is emitted by material located on a constant velocity 
surface, which will have a large geometrical extension (see e. g. 
Mihalas^ 1978, Ch. 14). However, the absorption of the O-star’s 
continuum can only arise in a narrow cross-section of the constant 
velocity surfaces, which lies along the line-of-sight between the 
observer and the 0-star. 

Figure 3 illustrates the geometrical setting of this 
discussion. The wind expands spherically about the WR core, which 
is at the center of the (p,z) - coordinate axes system. The degree 



Pi4&« 13 

OF POOR QUALITY 



Observer 


FIGURE 3; Simplified interpretation of the geometric setting giving 
rise to the phase dependent variations. 
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to which the absorption occurs at each frequency depends upon 
the optical depth along the line-of-sight to the 0-star at, that 
frequency, and this depends entirely upon the opacity of the material, 
y.(r) , and the velocity gradient in the wind. That is, 

T . K-M . XlE) 1 

dVz u i f dVCr) \ + V 

dz \ dr / t 


where V " V(r) (z/r), p - z/r, and r ■ (p^ + z^)^. The 
z 

opacity contains known atomic constants and the level populations: 


x(p) 


9 

ii e- 


mc 


(gf^) 
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where N^, are the poulacions of the lower and upper levels of the 

transition, respectively. The relevant velocity gradient in equation 

1 is the change in the 2 -component of the expansion velocity in 

the 2 direction (i.e., along the line of sight). 

Consider an orbital phase such as that illustrated in 

Figure 3. If we assume that the emitting material extends out to 

a distance z * Z fflax , all the material along the p*P * constant 

ray, from z * 0 to z * +Z ., will absorb radiation from the O-star 

max 

at wavelengths ranging from A 0 to A Q ( 1 - ^ V Z (P, Z’ )). Hence, 

A Q is the wavelength of the transition in question in the rest 
frame of the WR star. Therefore, this will generate a shortward- 


33 


shifted absorption. However, the material having a z-component 
velocity toward the O-star (i.e., V 2 (-z)) can also absorb, and will 
do so at wavelength longward of X Q . The maximum "redshifted" 
absorption will have a wavelength X Q (1 ~~V Z (P, -Zu)) where -Zc 
corresponds to the closest point to the O-star at which the WR 
wind structure is unmodified by the companion's presence. 

As the orbital phases varies, the impact parameter P does 
too, which leads to a vf-rying range in V z values. Thus, as one 
proceeds from orbital phase 0.75 to 0.0, wind material with a 
growing rang© of V values occults the O-star, and thus a progres- 
sively broader absorption is predicted. 

We must note that because we are dealing with ionic 
transitions, the opacity does not depend solely on the overall 
density dis tribution in the wind but also on the excitation conditions 
Even if the ded^ify and temperature distributions decrease mono- 
tonically with distance, the population at the lower level of the 
transition (which gives rise to the absorption) may not. Hence, the 
opacity for the given line may have a maximum at a certain distance 
from the stellar core, and this distance may vary for different 
transitions, depending on the excitation potentials. In Chapter III 
we will derive a non-mono tonic optical depth distribution. 

Alternatives to the selective atmospheric eclipses inter- 
pretation have been proposed, and rely on the presence of streams of 
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material or on a ’’heating” effect in the V?R hemisphere facing the 
0-star. In the first scenario, gas streams flowing through the 
Lagrangian points of the gravitational equipotentials of the system 
(Sahade 1958; Kuhi 1968) produce enhanced absorption at the 
appropriate orbital phases. However, WR + OB binary systems do not 
comform to a simple Roche Lobe, overflow geometry. First of all, 
the large wind velocities prevent the equipotential surfaces from 
having a significant effect on the mass outflow (Castor 1970b). 
Secondly, we now know that all OB stars have considerable mass out- 
flows of their own, leading to a wind-wind collision in the region 
between the WR and OB members of binary systems. The energy re- 
leased in this collision (Prilutskii and Usov 1976; Cooke et al. 

1978) renders any interpretation in terms of "cool" streams of 
gas extremely unlikely. The second interpretation is Somewhat 
related to the possible wind-wind collision « Many of the variations 
observed in WR binaries have led the observers to infer the presence 
of "hot" material located between the two members of the binary 
system (Kuhi 1968; Khaliullin and Cherapashchuk 1975). We will show 
that there is indeed evidence for a wind-wind collision in our 
observations , but that the effects of this interaction are secondary 
to those of the atmospheric eclipses. 

Figures 4(a) and 4(b) show the phase dependence of the 
emission and absorption equivalent widths, respectively. V444 Cyg, 


FIGURE 4: Equivalent widths of the (a) emission and (b) absorption 

components of N wV 1718 as a function of orbital phase. 

In (c) we plot the ratio of absorption to emission equiva- 
lent width as a function of phase, for comparison with 
this ratio in HD 193077, represented by the filled-in 
square. Other symbols correspond to close WN + OB systems 
and represent: . V444 Cyg, o ED 211853, x HD 90657 and 

+ HD 186943. 
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S3) 90657, HD 211853 and HD 186943 are repesented on the same scale, 
with different symbols. The similarity in the amplitude is 
striking, and apparently reflects the similarity of the WN winds 
in these cases. The gradual decrease in as the 0-star come 
out from behind the WR wind, accompanied by a gradual increase in 
emission indicates that the variations observed in the emission 
component are directly related to those in the absorption component. 

The ratio as a function of orbital phase is shown in 

Figure 4(c). The filled-in square represents thie value of this ratio 
derived from the HD 193077 spectrum. Since perturbations due to a 
close OB companion are minimal in this source, its W^/Wg value is 
closest to that expected from an unperturbed WR P Cygni feature, 
and can be used for comparison with the other sources. Thus, we may 
conclude that in the four variable sources, the N IV 1718 line is 
altered by enhanced absorption during orbital phases 0.8 < cp < 1.2, 
and perhaps reduced absorption at the opposite phases (i.e., <p ~ 0.5). 
The ratios are equal to the HD 193077 value at elongations, which 
tells us that radiative interaction effects at these phases are 
negligible. 

The variations resulting from selective atmospheric eclipses 
can be modeled using the computer code described in the Appendix, 
and compared with the observational data. In Figure 5 we have 
constructed an £ - a plot for the N IV 1718 line using the data 
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Comparison of model derived and observed (a, £) values. 

Curves marked Ml, M2 and M5 correspond to models with L ■? L 
'intensity ratios of 1:1, 1:2, and 1:5, respectively. ** 
Different symbols represent the different WN targets: 

. V444 Cyg, o HD 211853, x HD 90657, + m 186943, * HD 94546, 
A HD 193077 and X HD 50896. A£ and Aa are the propagated 
ivcertainties in the £ and a measurements, respectively. 
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in Table 2. As in the Appendix, e ■ F E /F C , and a * F^/F^, w ^ ere 
F c is the measured continuum level in the region of the emission 
line. Different symbols represent the different sources. For a 
given WIT + OB binary system, there is a large spread in a values, 
which results from the phase-dependent variations of the absorption 
component. Thus, at $ ~ 0, a* 0.4 - 0.5, increasing to a * 0.7 - 0.8 
at elongations and approaching unity (i.e., F^ 55 F^) at <J> ** 0.5 - 
0.6. The encircled "X" symbol on this plot represents HD 50896 
(SWP 13844 + SWP 14136) , whose value differs considerably from those 
of the other systems (including HD 193077) . HD 50896 is a system 
containing a WN5 Wolf-Rayet and a degenerate companion thought 
to be a neutron star (Firmani et al. 1980). In this case, the 
continuum level can be regarded as arising in the WN5 component of 
the system. In addition, the N IV 1718 line does not undergo 
variations like those we have described for the WN + OB systems 
(Rumpl, private communation) . 

We have used the (a,e) position of HD 50896 to determine the 
value of opacity at the surface of the core (CO) and the relative 
intemitiy of the continuum with respect to the emission line (PRIM) 
according to the procedure outlined in the Appendix. The computer 
code was run with these quantities, the additional parameters as 
listed in the Appendix (i.e., constant velocity law, etc.) and 
WN : 0B intensity ratios of 1:1, 1:2, and 1:5. 
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The curves in Figure 5 corresponds to the resulting model 
values of (a»e) pairs for orbital phases ranging from just past 
elongation (4> » 0.75) to just before geometrical eclipse of the 
0-star by the WR (<J> * 0.0), and are labelled according to the 
relative intensity of the 0-star. Near elongations, a -0.7 - 0.8, 
while near <j) * 0, a * 0.4, in very good agreement with the data. 

The dependence of £ on the relative intensity of the 0-star is 
evident, and the best agreement with the data is achieved with the 
M5 model. This corresponds to a 1:5 intensity ratio, as expected 
from optical eclipse observations of 7444 Cyg (Kuhi 1963). However, 
model e values become rather insensitive to relative values ’of 
WN : OB greater than 1:5, and ‘thus we can only set a lower limit 
to the intensity ratio. 

In conclusion, we show that by a) adopting wind parameters from 
HD 50896, b) including a luminous companion, and c) assuming selective 
atmospheric eclipses of the companion by the WN wind we can repro- 
duce observed (a,£) values for orbital phases 0.8 < <J> < 1.2. We 
are unable to generate the values of a observed at phases cp - 0.5, 
when our model does not take into account any "heating" effects or 
the finite diameters of the two opaque, continuum-emitting cores. 

In addition, we have assumed that the spectral classification 
represents actual physical conditions in the winds, and that there- 
fore we can apply informaion derived from one system (i.e., HD 50896) 
to other, slightly different systems. 
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In summary, we have presented results of empirical and semi- 
empirical tests, both leading to the conclusion that selective 
atmospheric eclipses provide a satisfactory explanation for the 
dominant phase dependent spectral variations observed in WR + OB 
binaries. We will now consider some of the second order effects of 
this interpretation, and in the next chapter obtain an optical 
depth distribution for the winds of our sample of WN stars. 

Wind - Wind Collisions 

In Figure 6(a) and 6(b) we have plotted the ratios of the 
images corresponding to "WR in front" divided by "O-star infront" for 
the five variable systems. Where possible, spectra corresponding 
to similar* phases have been co-added to increase the signal-to- 
noise factor. For HD 94546 it has been necessary to construct the 
ratio with the spectrum corresponding to $ - 0.5 in the numerator; 
i.e. , inverted with respect to the other cases . We note, however, 
that an error of 0.006 day in the orbital period, which is con- 
sistent with uncertainties in the period determination, could 
generate a phase shift of 0.5 when extrapolating zero phase from 
the initial epoch listed in Table 1 to the observations reported here. 

The changes in the spectra between phases near 0.0 and 0.5 
are evident in Figure 6. For example, the enhanced absorptions at 
4» — 0 show up here as depressions with respect to the level of 



FIGURE 6(a): Ratios of images corresponding to "WR in front" divided 

by 0— star in front". The images used in each case are: 
(SWP 15602 + SWP 15644) /(SWP 1558Q + SWP 15615) for V444 
Cyg, (SWP 15617) / (SWP 15629) for HD 211853, (SWP 15598)/ 
(SWP 15646) for HD 90657. Ratios are displaced vertically 
with respect to each other by 1.5 units. 
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FIGURE 6(b)*. Ratios of images as in Figure 6(a), corresponding to: 

(SWP 15578) /(SWP 15616) for HD 186943, and (SWP 15599)/ 
(SWP 15631 + SWP 15577) for HD 94546. 





unity. That is, if there were no variations, the ratios 

would be equal to unity at all wavelengths. Particularly striking 

0 

is the spike at -1540 A in Figure 6(a), which results from the 
narrower 3 IV 1550 absorption component at <f> - 0, as compared to 
(f> ~ 0.5. Since this line corresponds to a resonance transition, it 
can produce absorption even at extremely low densities (c. f. inter- 
stellar medium), and thus, a strong contribution from the O-star? s 
wind to the overall profile is expected. The fact that the 
absorption is broader when the 0-star is in front of the tTR implies 
that the 0-star's contribution occurs at higher expansion velocities 
than the contribution arising in the WR wind. This means that the 
0-star companions in the systems of Figure 6(a) have C IV out to 
higher wind velocities than do the WR's. Published results for 
terminal velocities in 0-stars (Garmany ef w , al. 1981) and in WR's 
(Barlow et al. 1981) suggest that the early 0-s tars do indeed 
have higher terminal velocities. Though adequate phase coverage is 
lacking, it may be significant that the spike is absent in the 
HD 186943 ratio, since this system contains the latest 0-type 
companion of our sample. 

Strong wind-wind interaction effects are suggested by the 
curious behavior of the shortwardmost edge of C IV 1550. That is, 
the fact that the high-velocity portion of the absorption is 
present primarily when the O-star is in front of the WR, and 
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diminishes when it is in back implies that the 0-s tar's wind 
structure is disturbed. In the case of a wind-wind collision, the 
interaction would affect primarily the outermost portions of the 
0-star* s wind, which is where maximum expansion velocities have 
already been reached. Because the outflow from the WR is more 
massive than that of the 0-star, the former is expected to 
dominate, and one can envision the O-star's wind being stopped, 
and the ionization structure of its outer portions altered by the 
collision. Thus, the spikes in Figure 6(a) provide observational 
evidence for wind-wind collisions in WR +* OB binaries . 

Figures 2(a) , (b) , and perhaps (c) suggest that the 
alteration of the O-star* s wind may not be limited to the hemis- 
phere facing the WR star. Xn these figures, the C IV 1550 absorption 
appears to be narrower close to elongation, as well as at <p - 0, 
though this may be a result of the relative orbital motions. 

Evidence from X-ray observations for the existence of wind- 
wind collision regions in WR binaries may exist (Sanders, Cassinelli 
and van der Hucht, 1982) since the X-ray spectral distribution of 
one source (y Vel) is harder (T - 10^ K) compared to other WR 
sources. The latter's spectral distribution appears to be con- 
sistent with the X-rays being generated by small shocks within the 
stellar wind (Lucy and White 1982). 



Possible " and Fe VI Contributions 


With the exception of V444 Cyg at orbital phases $ * 0.00 

0 

+0.03, no appreciable changes in the continua at X > 1350 A occur 

in our SWP images nor in LWR (Long Wavelength Redundant, XX 1900 - 

3300 A) spectra of HD 211853 and HD. 18943 (Hutchings and Massey 1983) . 

0 

However, changes shortward of 1500 A are significant, especially in 

systems of Figure 6(a). These changes are approximately of the 

same magnitude as those in N IV 1718 , and have the same phase 

dependence. Thus, it is very probable that they result from line 

variablity. That is, closely spaced, numerous lines all behaving in 

the same manner would, given their Doppler broadening, produce 

the observed effects. In fact, there is a large concentration of 

0 

lines at X < 1550 A, most of them due to transitions of Fe V and 

F® VI ions, as a search through published line lists (Reader et al. 

* 

1980; Ekberg 1975a, 1975b; Ekberg and Edlen 1978) can easily show. 
Specifically, Fe V has 236 lines within the wavelength interval 

o a 

XX 1300 - 1500 A, but only 78 weak lines between 1500 - 1700 A. The 

O 0 

strongest lines all lie between 1350 A and 1486 A. As far as Fe VI 

is concerned, it has no. transitions at wavelengths greater than 
0 

1500 A, with all the strong ones having wavelengths shortward of 
0 

1380 A. It is also worth noting the Fe V has the same ionization 
potential as N IV, thus making its existence in the WN4 - 6 winds 
very reasonable. 


47 


> 


In Table 4 we have listed approximate wavelength positions 
of all emission features visible on the calcomp plots in our 
observations along with suggested identifications. With respect 
to the possible presence of Fe V and Fe IV, we can only state that at 
every wavelength position corresponding to strong Fe V lines , 
there is an emission feature in our spectra. The low resolution is 
not altogether the limiting factor for appropriate line identi- 
fications, since the large intrinsic width of WR emission lines 
results in severe line blending. As shown from high-dispersion 
optical spectra (Bappu 1973), several identifications are usually 
possible for a given line. In our case, Si III and Si IV have 
emission lines coincidentwith most of the Fe V and Fe VI lines. 

0 0 o 

The exceptions are the lines at 1320 A, 1330 A, and 1586 A, where 
Fe V and Fe VI appear to be the only reasonable contributors, 
according to the line lists we searched. 

Lines of Fe V have been identified in the spectra of the 
binary UW CMa (Drechsel et al. 1981), two WR's (Fitzpatrick 1982), 
and an 0-type subdwarf (Bruhweiler 1981) , although in the latter 
three objects they are presumably photospheric absorptions. 

In conclusion, we attribute the apparent continuum variations 

o 

shortward of 1550 A in our IUE images to line variations, as 
observed in N IV 1718, with the major contributors being Fe V and 
possibility Fe VI. All lines having adequately populated lower 
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TABLE 4- Suggested Line IdeaEificationa 


* (A) 


1243 




N V 1243 

1287 


Si IV 1286 


Fa VI 1285-88 

1311 

Si III 1312 



Fa VI 1308,10 

1320 



Fa V 1318-23 


1330 



Fa V 1330 

Fa VI 1329-31 

1340 

Si III 1341-43 


Fa V 1344 

Fa VI 1338 

1364 

Si III 1362-65 


Fa V 1361-66 

Fe VI 1362 

1377 

Si III 1373 

0 V 1371 

Fa V 1373-76 

Fa IV 1371-75 

1395 


Si IV 1394 

Fa V 139* 


1407 


0 IV 1407 

Fa V 1407-09 


1420 

Si III 1417 


Fa V 1418-2* 


1431 

Si III 1433 

• 

Fa V 1430 


1440 

Si III 1442 


Fa V 1440-42 


1445 

Si III 1447 


Fa V 1447-49 


1470 

0 III 1477 

Fa IV 1473 

Fa V 1465-79 


1486 


N IV 1486 



1505 

Si III 1500-06 



0 V 1507 

1532 


Fe IV 1530-34 



155-3 


C IV 1548-51 


N V 1549 

1575 


Fa IV 1577-78 



1586 




Fe IV 1583 

1625 


Si IV 1624 

Fe V 1623-27 

N V 1620 

1682 


H rV 1687 



1705 


Fa IV 1705 


0 V 1708 

1721 


N IV 1718 



1732 

N III 1730 





49 


levels can partake in the absorption of the O-star' s continuum at 
relevant orbital phases, but only Fe V and Fe VI have the sufficient 
number of lines to generate a pseudo-continuum. 

Radial Velocity Curves 

One consequence of the selective atmospheric eclipses 
is the alteration of radial velocity curves. In Figure 7 we have 
plotted as a function of orbital phase the positions of maximum 
emission (Vg) and deepest absorption (V^) for V444 Cyg, HD 90657, 

HD 136943, and HD 211853. The mean velocity for each system has 
been subtracted in each case. 

In both 7(a) and 7(b) , minimum velocity occurs at <p - 0.75, 
as it should for a true radial velocity curve, since this is when 
the WR is approaching the observer and its lines have the maximum 
shortward shift. At $ - 0.25, the velocity should be maximum. 
However, although the absorption component in Figure 7(b) does 
describe a somewhat typical RV curve, the emission component 
definitely does not. In fact, it has a maximum at <J> - 0, in 
addition to its amplitude being larger by about 50^. This is an 
example of the well-known problems besieging RV curves in WR stars; 
namely, for a given binary system, the RV curves derived from 
different emission lines have phase shifts, they have different 
amplitudes and zero-velocity values, and sometimes are distorted so 



<V*> iio* km V 4 -<V 4 >U0 2 km*-') 



iraa (b) 
ital ph 
Figure 
ties. 
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as ta suggest non-circular orbits. This may now be understood in 
terms of the atmospheric eclipse interpretation. 

At <p - 0, the WR wind is absorbing the 0-star’ s radiation. 

Since the latter. is behind the MR wind, absorption will occur both 
shortward and longward of line center. However, due to wind-wind 
interaction effects the strongest absorption will occur shortward 
of line, center., This will drive the centroid of the emission 
component to longer wavelengths at <p - 0, as we see in Figure 7, 
distorting it. At the opposite phase, the absence of the shortward- 
shifted absorption will restore the emission line to its rest wave- 
length. The net result is a phase shift with respect to the true 
RV curve. The degree to which’ this effect is important depends upon 
the orbital inclination of the system, its separation, and the 

t 

relative luminosity of the OB component. With respect to a given ? 

system, lines which show phase-dependent profile variations will 

most likely produce modified RV curves. Thus, for the determination ; 

i. 

it 

of orbital parameters, it is best to choose emission lines that do ? 

not generally have a shortward-shifted absorption, whose zero- 
velocity best matches that of the OB component, and whose RV curve 
best resembles a true RV curve. These are the criteria generally 
applied, when possible. 
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‘ Chapter III 

EMPIRICAL OPTICAL DEPTH DISTRIBUTION 

la Chapter I we pointed out that the observed emission line 
profiles In WR stars can be fit by model profiles with dissimilar 
combinations of source function and optical depth distributions 
(Castor 1970a). Thus, the wind structures cannot be unambiguously 
derived through theoretical model fitting alone. Therefore, an 
empirical determination of the optical depth distribution would be 
very useful in constraining the free parameters in the model 
calculations. 

The optical depth along the line-of-sight to the 0-star at 
the wavelength corresponding to the minimum in the N IV 1718 
absorption (i.e., V^) can be defined as 



where F^ (<f>) is the measured flux at orbital phases 0.75 < <f> < 1.25 
and Fq (c}>) is the corresponding measured continuum, flux* The factor 
0.2 is substracted to account for the contribution- to the absorption 
at V A resulting from the projection of the MR wind onto the WR core. 
This value is taken from measurements of F a /Fq in HD 193077 and HD 
94546 and at orbital phases -0.7 in HD 186943, which is the widest 
of our close WN + OB binary systems. 
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We will derive the dependence of optical depth upon impact 
parameter, p, which is related to the orbital phase by 

p * a (sin^ (2 tt<{>) + cos^ (27r<|)) cos^ i]* 1 

where a is the orbital separation, and i is the inclination of 
the orbital plane. 

Figure 8(a) is a plot of as a function of impact parameter 
for values listed in Table 5. Different symbols correspond to 
different binaries, as in earlier figures. Note that by combining 
results from different systems in this manner we are assuming that 
the wind structures in WN4-6 stars are basically the same. The 
curve shown in this figure corresponds to a least-mean-squares fit 
to the data, resulting in 



Before drawing any conclusions, the following cautionary 
remarks may be made. First, the material responsible for the 
absorption at V A has a net outward (toward the observer) velocity. 
Thus, it has coordinates (p,z y ) with z > 0 (see Figure 3). We 
have no means of determining z. Therefore, the dependence of T 2 on 
p can only be regarded as an approximation to the true t(r) 


FIGURE 8: (a) Optical depth distribution as measured from the minimum 

flux level in the N IV 1718 absorption component as a 
function of impact parameter. A least-mean square fit to 
the data at p > 14 R fi is illustrated. 

(b) Optical depth distribution as measured at X 0 , the 
rest wavelength of the N IV 1718 transition in the WR 
frame. A power law of the form t 0 «p"^ is superposed. 
Propagated uncertainties for l Q values are At 0 * ±0.4 
(see page 6S . 
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TABLE 5: Impact Parameters and Computed Optical Depths. 


SWP 

PHASE 

p<V 

T 2 

T 0 

V444 Cyg 





13832 

0.917 

19.0 

0.51 

0.95 

15380 

0.852 

29.0 

0.17 

0.48 

15387 

0.992 

7.6 

0.77 

0.71 

15388 

0.998 

7.4 

0.72 

0.60 

15389 

0.003 

7.5 

0.77 

0.51 

15390 

0.008 

7.6 

0.77 

0.72 

15391 

0.015 

. 8.1 

0.62 

0.74 

15392 

0.020 

3.6 

0.72 

0.82 

15394 

0.030 

9.9 

0.62 

0.56 

15396 

0.051 

13.5 

0.77 

0.34 

153S7 

. 0.056 

14.2 

0.64 

0.87 

15399 

0.066 

15.9 

0.64 

0.94 

15400 

0.071 

16.8 

0.60 

1.01 

15401 

0.076 

17.7 

0.62 

0.95 

15402 

0.081 

18.6 

0.53 

0.86 

15403 

0.085 

19.3 

0.64 • 

0.78 

15406 

0.103 

22.4 

0.56 

0.83 

15407 

0.108 

23.2 

0.53 

0.65 

15408 

" 0.114 

24.1 

0.47 

0.73 
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TABLE 5; Impact Parameters* and Computed Optical Depths (continued) 


SWP 

PHASE 

p<V 

X 2 

T 

0 

15588 

0.835 

31.0 

0.29 

0.72 

15602 

0.073 

'17.2 

0.53 

0.88 

15644 

0.974 

9.4 

0.49. 

0.84 

HD 90657 

15584 

0.83 

49 

0.19 

0.23 

15S98 

0.95 

26 

0.40 

0.54 

15618 

0.19 

. 51 

0,11 

0,0 


HD 186943 


15578 

0.11 

47 

15586 

0.22 

67 


G.Z 


0.0 


0.2 
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distribution* The accuracy of the approximation depends upon z 
being small enough to justify writing r 2 * p 2 + z2 * p2, Second, at 
small impact parameters the finite dimensions of the stellar cores 
become important, and they have been neglected. Third, the 
assumption of similarity between winds in WN4-6 stars may be 
questionable, and the points at p > 30 R^ correspond to the two 
WN4 stars, while those a p < 30 R^ belong to the WN5. 

An opti cal depth distribution of the type t « r is 
consistent with a constant wind velocity distribution. That is, 

for a constant vfr) - V w , conservation of mass flow in (M « r 2 pV) 

-2 

demands p « r , which if we assume that the populations of the 
levels go as the total density, implies that 



where we have used equations .1 and 2 with the assumption of z < p. 

Thus, with previously stated reservations, we conclude that, to a 

first approximation, the velocity distribution in the winds of 

early WN stars is constant, at least beweeen -14 and 66 R . 

Cherepashchuk, Eaton and Khaliullin (1981) concluded that the 

velocity is constant at r > 12 R from their analysis of the 

a 

optical eclipse light curves. Thus, the consistency of our results 
with previous observations is encouraging. 
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It is interesting; to note that i£ > 0 for impact parameters as 
large as 66 R fl , which means that the N IV ionization stage exists 
at least out to this distance from the WR core. Castor and van 
Blerkom (1970) found that the He II emitting region" in the WN6 
star HD 192163 extends out to 70 R^. Assuming that this result is 
applicable to our WN4-6 stars, we are led to conclude that ions 
with different ionization potentials do indeed coexist in WN winds, 
as concluded by Willis (1982) . 

We will now describe an intuitively attractive method by 
which the measured optical depth can be associated with an exact 
position in the wind, thereby deriving the true optical depth 
distribution as a function of radius. 

Consider the rest wavelength (1 0 ) of a P Cyg feature such 
as N IV 1718, and consider an orbital phase 0.8 < <{> < 1.2; i.e., 
when the 0-star ia behind a portion of the WR wind (z < 0 in Figure 
3). Regardless of the velocity distribution in the WR wind, the 
only contribution to the absorption at the rest wavelength will come 
from material flowing perpendicular to the line-of-sight to the 
0-star. This material is- situated at a distance r * p from the 
WR core, and at 2 * 0 + ** 2 th* wl ^ ere z th is t ^ e tbermal length; 
i.e., the distance over which V 2 changes by v^, the thermal 
velocity of the particles at that point. V_ is the z-component 
of the expansion velocity. That is, a photon with wavelength \ Q 
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can be absorbed only by material traveling with a relative 
velocity (with i,“«pect to the source) v < v^. Since in an 
expanding atmosphere t q counts only the absorbers in a narrow 
cross-section of the column along the line^a£ sight (according to 
the Sobolev Theory) , it can be expressed as 


/th/2 

T 0 - Icp^rMz 

" Z th/2 


where p, is the effective density of absorbing ions, and K is the 
i. 

cross section for the transition. 

A thermal, path length z til can be derived in the following 
manner : 


V 2 (r) - v(r) 


z 

(p2 + z ^) 1 


v 


th 


Z “ z 


th 


since V,_(p,0) • 0. Solving for z we find 


P V th 

th * 


1 


z * z 


v(r) 


C 1 ' v th/v Z (r)) % 
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For typical WE winds, v(r) » iOO Km and v^ < 30 Km s“^ 
(T e < 50,000 K) . Therefore, ( v t t/ v ^) 2 <<: and 


Z Lh 


P v 


th 


v(r) 


P v th 
v(p) 


since 2*0, and rap. 

Hence, the optical depth at the rest wavelength of an 
emission line can be expressed as 


Kp^Cp) 


PVth 

v(p) 


3 


where 


Kp (p) ■ X(p) 
i 


2 

TTe 


me 





(Hihalas, p. 479). Equation 3 could also have been obtained directly 
from eqs. 1 and 2 by setting z * 0. 
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If the density of absorbing ions followed the total density 
distribution, f^(p) would be constant, and the above equation 
could be used directly to derive v(r) by measuring t q * The impact 
parameter p is presumably known. However, we cannot guarantee that 
f i (p) * constant, which would mean that the ionization and 
excitation levels are frozen in the wind. As mentioned previously, 
there is evidence for a frozen ionization structure in WN winds, 
but not for a frozen excitation structure (Willis 1982) . 
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For illustrative purposes, we will now proceed to *pply the 
method we have just described to our IUE observation, thou- \ any 
results that may be derived can only be regarded as tentative, 
due to the low resolution of the data. 

We wish to determine the opacity of X q due to the bound 
bound transition 2S2P 1 ? - 2p 2 in N IV resulting from the 
absorption by ions in the WN wind of the O-star' s continuum. By 
definition, radiation which emerges from an absorbing medium is 
expressed as 


F * F e" T 
o 

where F 0 is the incident radiation. In this case, F Q corresponds 
to the continuum of the O-star. Given the errors which will be 
introduced by the low resolution of our observations, we are 
justified in asstuning that the O-star continuum dominates the UV 
continuous energy distribution and that the measurable continuum 
level is generated primarily by the O-star. This is not a bad 
assumption since optical light curves of V444 Cyg suggest that 
the UV luminosity ratio is 5:1, with the O-star being the brighter. 
Hence, we adopt F^ » F , where F, is the measurable continuum level. 

r Q C C 


A s long as the 0«star is far enough away from the material 
flowing along the p axis, the interaction effects should not perturb 
the material emitting at Thus, we assume that the emission 

at X Q is constant throughout the orbital phases we will consider. 
That is, we assume that the WR wind is stable, and that all alter- 
ations of the profiles are due to underlying absorption. 

At orbital phases <p - 0<,5, the flux we measure at X Q can be 
approximated by 


F , 0,5) - F (0.5) + P ft ) 7. 

mo c e o 

since the interaction effects should alter the profile mostly 
in the shortwardmost regions. When the 0-star is behind a part 
of the wind at phase (p, the flux we measure is 

F m C\ 0 , « - F c «) e" T ° + F e a o ) 8 

-T 

where F c (<£) =* F c (0.5) e c and^r^ takes into account any 
continuum opacity sources (such as electron scattering) , We will 
take F c (<$») tot be: the measured, continuum level at orbital phase (p, 

o' 

at the 1718 A region. 


Combining equations 7 and 8 and solving for t q we obtain 


'In 


fjk±°L 


<ft) + F e (0.5) - F 
F c (« 


- <* n - °~ 5 > 1 


The values of F m were interpolated from a least mean squares 
fit to five points* on * the., profile, centered on the rest wavelength 
of N IV 1718. S in the WR ? s reference frame. The current IUE data 

o 

extraction, procedures permit flux values to be obtained at 1.2 A 
intervals. Thus, 5 points represent one resolution element. All 
profiles were inspected visually to verify the interpolation. 

The F c values were obtained by visual interpolation across the 
? Cyg feature, with the RDAF routine FEATURES. 

Figure 8(b) is a plot of T 0 as a function of impact para- 
meter for the same spectra represented in Figure 8(a). The 
propagated uncertainties in this case are very large (At q ~ 0.4) 
and are not plotted. Despite this, we may venture a few comments. 
First, there is again a decrease in optical depth for r < 14 R^ 
as in Figure 8(a), while the data at larger radii can be represented 
by a power law of the form T 0 * r Second, the rapid decrease 

in optical depth at small distances from the core is reminiscent 
of the optical depth distributions employed by Castor (1970a) to 
produce realistic model profiles. Castor suggests that different 
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excitation conditions may prevail close to the MR core, as appears 
to be the case according to Figures 8(a) and 8(b). 

Equation 6 can be used to study the excitation structure of 
the wind, once the velocity distribution is derived. When investi- 
gating a given binary system, the ratio 

f i (p) V 2 (Pf) ^ T <p) P 

fiCPf) V 2 (p) T(p ? ) p f 

can be constructed, thus eliminating the constanta, anxi obtaining the; 
density of absorbing ions relative to that at p f . In our case, we 
have fixed p^ * 31 R^ and plotted the points corresponding to 
f‘ i (p)/^ (31) in Figure 9. Me are adopting the result of the T 2 
analysis; i.e., v(r) * constant and using only V444 Cyg data. This 
figure seems to imply that the population of the. lower level of the 
N IV 1718 transition at r * 7 - 10 R is a factor of two smaller 
than at r - 14 R Q . Thus we may speculate that the population of the 
lower level of the N IV 1718 transition increases with distance from 
the MR core. 

Willis (1982) has established a correlation: between the 
excitation potential (EP) of the lower level of a transition and 
t;he velocity of the deepest portion of the. absorption. (V^, in our 
notation). He finds thac for a given star, increases for lines 


«r)/&) 
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6. IQ II. 18. 22. 26. 30. 


r/Ro 

FIGURE 9: Possible dependence of the population of the N IV 1718 

lower level on distance from the WR core, suggesting a 
changing excitation structure as a function of distance. 


mtecr. 
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with lower EF* Thus, if one assumes an accelerating flow, 
excitation conditions decrease in the wind with increasing distance. 
According to our results, and those of Cherepashchuk, Eaton and 
KhaliuLXin. (198Q) » the accelerating region must lie close to the 
WR core. The increasing optical depth between r - 7 and r - 16 R^ 
in Figure 8(b) is in agreement with Willis’ velocity - EP 
correlation since the EP of the N IV 1718 lower level is small 
(16.5 eV)'. That is, decreasing excitation conditions in the wind 
translate into increasing populations of low-lying levels. 

We conclude by stressing that, due to the limitations of our 
data, the last remarks may only be regarded as speculative at this 
time. 
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Chapter IV 

CONCLUSIONS AND RECOMMENDATIONS 

The objectives of this investigation have been to analyze the vari- 
ability of Wolf-Rayet binaries (WN + OB) in the ultraviolet spectral 

0 

region XX 1200 -1900 A, and determine whether the periodic, phase- 
dependent variations can be used to empirically derive the physical 
properties of the WN winds. We have found* the- variations in the UV to 
be similar in nature to those observed in the optical spectra. That 
is, there is a strengthening of absorption: components in P Cygni- 
type features at orbital phases in which the 0-star is behind the WR 
wind. With the aid of a computer code which models this type of 
variations, the dominant mechanism producing the variations is 
shown to be selective atmospheric eclipses of the OB component by 
the WR wind. Based on this explanation^ a straightforward technique 
is proposed as a first step in determing the wind structure. The 
analysis of IUE observations applying this technique results in an 
empirical determination of an optical depth distribution in the WN . 
winds which implies that a constant velocity law prevails at distances 
of 16 < r <66 R^. In addition, it suggest that the N IV 
ionization stage exist at least out to 66 R consistent with the 

9 

proposed constant ionization structure (Willis 1982). However, there 
may be an indication of a varying excitation structure at r < 16 R # , 
which would be consistent with the observed velocity-excitation 
potential correlation which has been derived (Willis 1982) . 
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The next steps that must be taken Include: a) modeling the 

entire P Cygni profile using empirically determined parameters, 
b) applying the analyses to emission lines of different ionization 
and excitation potentials, and c) obtaining full phase coverage of 
additional sources. 

As Castor (1970a) pointed out, WR profiles can be modeled and 
reproduced with a variety of dissimilar combinations of optical 
depth and source functions. We propose that the additional infor- 
mation required to restrict the free parameters in model profiles may 
be obtained from optical depth analyses similar to the ones we have 
performed. In particular, optical depth distributions such as those 
in Figure 7 can be used to fit model profiles to the observations. 

Evidence supporting the existence of wind-wind collisions is 
provided by the variations in the C IV resonance line absorption 
component. Thus, an additional step that must be taken from the 
modeling standpoint is to include the wind-wind effects into the 
code. It is not clear at this point how one is to proceed on 
this question, unless the general solution to the problem of wind- 
wind collisions is derived. Furthermore, the interaction not only 
affects the OB star's wind; the structure of the WR wind must also 
be modified to a certain extent, and hence, a non-spherically 
symmetric wind velocity distribution would be more realistic for 
actual model profile fitting attempts. 


Finally, it is important to repeat the analysis of Chapter III 
on additional emission lines having different ionization and excit- 
tation potentials. Excellent candidates are the lines of N V 4603, 
4619, and N IV 3483. These undergo periodic variations of the same 
nature as N IV 1718 and may provide the additional information needed 
to separate the excitation and velocity structures. 
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APPENDIX 

A computer code was constructed to model the phase-dependent 
profile variations in binary systems in which one member has a 
spherically expanding wind, or envelope. No interaction effects 
other than the absorption of radiation are taken into account. 

The treatment of the radiative transfer is based on the Sobolev 
Theory (1958, 1960) . Input parameters include the velocity and 
opacity distributions, which are assumed to be monotonic functions 
of distance. The source function, is taken to be constant. 

The profiles we will consider are composed of three parts. 

a) a pure emission line arising in the spherically expanding 
wind, disregarding the presence of a continuum emitting 
core. 

b) a shortward-shif ted absorption component produced by 
the absorption in the wind of the continuum radiation 
emitted by the core (WR)* 

c) an absorption of the companion's continuum radiation by 
material in the WR wind at emission— line frequencies. 

The emission (a) with the absorption (b) constitute a standard 
P Cygni profile, where the absorption component ha3 a net shift to 
shorter wavelengths. This absorption arises in the column of wind 
material which is projected upon the central continuum-emitting core, 
and lies along the line of sight to the observer. Iir Figure 3, this 
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material is flowing along the +z axis, at p « 0. Clearly, for 
an expanding velocity field, the absorption can only occur at wave- 
lengths shorter than X q , the wavelength of the transition in the 
rest frame of the central core. The depth of the absorption will 
depend upon the combination of opacity, velocity gradient along 
the line of sight, and intensity of the continuum source. 

The only difference between the absorptions in (b) and (c) 
is that the latter can occur at X ^ since the secondary continuum 
source is behind a column of material expanding both away from it 
(i.e., toward the observer) and towards it. 

It is. worth noting that given our assumptions, the absorptions 
affect the emission line only insofar as they lower the underlying 
continuum level at. selected frequencies. That is, the emission line 
retains its intrinsic*" profile and intensity, but the superposition 
of the absorptions* produces an apparent decrease in emission 
intensity. 

The free parameters in the model are: 

RMAX maximum extent of the emitting region. This is a 

useful parameter since it permits a cut-off in the 
emitting region before the density has fallen off 
considerably, which allows the possibility of a 
changing ionization structure. 
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VO,VN specify the power law velocity distribution v(r) * 

VO t W 

CQ,CN specify the power law opacity distribution CHI(r) * 

00 r® 

PRIM, SEC relative intensities of the WR and 0-star's continua 
RAD, VEL orbital radius and veloctiy, respectively 
WIDTH total width of SEC's emission. 

All distances are defined relative to the smallest radius of the 
WR emitting region, which may be taken as the WR core radius, R c . 
Velocities are given in thermal units, and for convenience we have 
chosen V^ * 10 km s”^. Intensities are in arbitrary units. 

The large number of free parameters can be reduced to a great 
extent by adopting the results of Chapter III and using known orbital 
parameters. Hence, we may fix. tha following^values : 

VO * 150. approximate terminal velocity 1500 km s“^ 

VN * 0 constant velocity 

CN * -2 assume P^(r) 11 P(r) 

VEL * 30. typical orbital velocity * 300 km s ^ 

WIDTH », 500. secondary's (O-star) emission is continuum 

. RMAX » 70/R c • maximum extent of emission from FIG 8 

RAD * 36 /R c orbital radius in V444 Cyg 

We have devised a method to reduce the free parameters further, 
through a parametrization of the P Cygni profiles and a comparison 
with easily measurable quantities. 
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We define the following quantities: 

e - V F c 

a - F a /F c 

where Fg, F^, F^ are equivalent to the measured quantities in Table 
2. If one chooses a velocity law, and specifies the inital opacity 
(CO), the code can be run for a grid of PRIM values. In the case 
of a single star, F c - PRIM, Fg » Fg + PRIM, and F^ » PRIM <e“ T ) . 

F £ is the maximum intensity of the pure computed emission. 

Using the resulting Fg, F^, and F c values, one can. obtain (a,e) 
points for every PRIM in the grid. When plotted (e vs a) , a smooth 
curve results. By choosing various values of CO, a family of curves 
is derived, as illustrated in Figure 10, where each curve is labeled 
with its corresponding CO value. It is Important to note that the 
curves do not cross, so that a given point on the plot specifies 
the value of CO. The encircled "X” represents the (a,e) point for 
the N IV 1718 line in HD 50896, and thus provides the value of CO 
to be used in model profile computations assuming, a constant velocity 
distribution. It also provides the value of PRIM, since PRIM =* 
Fg/(e-l). Fg is determined by the code once the opacity and velocity 
distributions are specified. 


FIGURE 10: Model (a,s) values for a single WE star as a function of 

WR continuum intensity (PRIM) and initial opacity at the 
surface of the core (CO) for a constant velocity law. 
Curves are identified according to the CO value. The en- 
circled "X" represents the N IV 1718 line in HD 50896. 
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